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unitary transformations: two density matrices are locally equivalent if and only if all these invariants 
have equal values. 
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Introduction. - As a fundamental phenomenon in 
quantum mechanics, the quantum nonlocality has been 
recently extensively investigated. Nonlocally quantum 
correlated states, like quantum entangled states [l| or 
states with nonzero quantum discord [2|, play very im- 
portant roles in many quantum information processing 
such as quantum computation [3| , teleportation [4J , dense 
coding 5], cryptography fgj and assisted optimal state 
discrimination [7]. 

Due to the fact that the nonlocal properties, e.g. the 
quantum entanglement of two parts of a quantum sys- 
tem, remain invariant under local unitary transforma- 
tions, they can be characterized in principle by the com- 
plete set of invariants under local unitary transforma- 
tions. For instance, the trace norms of realigned or par- 
tial transposed density matrices in entanglement measure 
[|| and the separability criteria [§] are some of these in- 
variants. A complete set of invariants gives rise to full 
classification of the quantum states under local unitary 
transformations. 

Many approaches to construct invariants of local uni- 
tary transformations have been pre sented in recent years. 



The method developed in [lOt in principle, allows 
one to compute all the invariants of local unitary trans- 
formations, though it is not easy to perform it opera- 
tionally. In 12j, a complete set of 18 polynomial invari- 



ants is presented for the locally unitary equivalence of 
two qubit mixed states. Partial results have been ob- 
tained for three qubits states 13j, some generic mixed 
states 14 - lj|] , tripartite pure and mixed states 17 1 . The 



local unitary equivalence problem for multipartite pure 
qubit states has been solved in 18]. Then the problem 
for arbitrary dimensional multipartite pure states are also 
solved recently [191 ] . However, for mixed sates, generally 
we still have no operational criteria to judge the equiva- 
lence of two arbitrary dimensional bipartite states under 
local unitary transformations. 

In this letter, we study the nonlocal properties of arbi- 
trary dimensional bipartite quantum systems and solve 
the local equivalence problem by presenting a complete 



set of invariants such that two density matrices are lo- 
cally equivalent if and only if all these invariants have 
the equal values in these density matrices. These in- 
variants can be explicitly calculated and give rise to an 
operational way to judge the local unitary equivalence 
for nondegenerate density matrices. For degenerate case, 
due to the eigenvector decompositions of a given state 
are not unique in the degenerate eigenvectors' subspace, 
the approach is no longer operational since the expres- 
sions of the set of invariants are not unique. Nevertheless 
the set of invariants is still complete in the sense that if, 
in a suitable eigenvector decompositions, two states have 
the same values of the invariants, they must be equiva- 
lent under local unitary transformations. In particular, 
we also present a set of invariants that are independent 
on the detailed eigenvector decompositions of a quantum 
state. These invariants present a necessary criterion on 
local unitary equivalence. 

The linear space formed by the invariants. — We first 
establish a linear space spanned by matrices whose traces 
are the invariants under local unitary transformations. 

Let H be an iV-dimensional complex Hilbert space, 
with i = 1, ...,N, an orthonormal basis. Let p be a 
density matrix defined on H®H with rank(p) = n < N 2 . 
p can be generally written as 



^2\\Vi)(Vi\, 



(1) 



»=i 



where \vi) is the eigenvector with respect to the nonzero 
eigenvalue A.;. is a normalized bipartite pure state of 
the form: 



iV 



h) = a ki\ kl )i a li g 
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where * stands for complex conjugation. 

Two density matrices p and p' are said to be equivalent 
under local unitary transformations if there exist unitary 
operators U\ (resp. U2) on the first (resp. second) space 



2 



of H (g> H such that 

^ = (Z7i ® I7 a )p(I7i <8> I7 2 )+, (2) 

where f denotes transpose and complex conjugation. 

To solve the local equivalence problem of the density 
matrices p and p' , it is sufficient to find the complete set 
of invariants under local unitary transformation Ui®U 2 - 
Lemma 1. The following quantities are invariants under 
local unitary transformations: 

J s {p)=T ri {Tr 2P % 
Tr[(^j)(^t)...(^4)], 

where s = 1,...,N , i,j,k,l,...,h,p = l,...,n, Tr\ and 
Tr2 stand for the traces over the first and second Hilbert 
spaces respectively, Ai denotes the matrix with entries 
given by {AA k i = a l kl , i = 1, 2, • • • , n. 
[Proof], Let U\ and U2 be unitary transformations sat- 
isfying U x Ul = U\Ui = U 2 U ] 2 = UlU 2 = 1. Un- 
der the local unitary transformation U\ ® U2, we have 
p — > p' = Ui ® U2 p U\®U\. Correspondingly, we have 
\ui) — > \uQ = Ux®U2\vi), or equivalently Ai is mapped to 
A- = U1A1U2, where t/J is the transpose of U?.. Therefore 

AiA't = U.AApl A?A>=u;AtA 3 Ut (4) 

for any i,j = l,...,n. By using relations (|4]) 
and taking into account that Tr2\vi)(vi\ ~ AiA], 
it is straightforward to verify that J s (p') — 
Tn[Y? i=1 KTr2{\v'Mm = Tn^XfA'A?} = 
J s (p), Tr[{A' i A'}){A' k A^)...{A' h A^)] 

Tr[U x {AiA*){A k 4) ■ ■ ■ (^4)^] 

Tr[{AiA^){AkA\) ■ ■ ■ (A h Ap)]. Hence the quantities in 
([3]) are invariants of local unitary transformations. □ 

If two density matrices are equivalent under local uni- 
tary transformations, then their corresponding invari- 
ants in ([3]) must have the same values. Before prov- 
ing that if two density matrices have the same val- 
ues of all the invariants in ([3]), they must be equiva- 
lent under local unitary transformations, we first claim 
that the set of independent invariants in ([3]) is finite. 
In fact, each factor A4A*, belongs to the finite dimen- 
sional algebra Mat(N), which has a linear basis Eij, 
i,j = l,--- ,N. Therefore, denoting r the number of 
factors like {A,A]) in Tr[(AiA})(A k A}) ■ ■ ■ {A h Al)\, we 
have that r is at most iV 2 . Subjecting to the variations 
of the subindices i, j, k, I, h,p, there could be many in- 
variants of this form. However, r may be much less than 
./V 2 for given i, j, k, I, h,p. For instance for the case of 
i = j — k — I... = h = p, T takes values from 2 to N 
(the case r = 1 is trivial since Tr(AiAj) = 1 due to the 
normalization of the state \vi)). The same is true also for 
the invariants Tr[{A\Aj){A\A{) ■ ■ ■ {A\A p )\. 



Denote the subalgebra of Mat(N) spanned by prod- 
ucts of AiA'j by 

M(p) = 8pan{(AiA$)(A k A]) ■ ■ ■ {A h A\)}, 

i,j,k,l,...,h,p = 1, ...,n. Obviously IR(p) is a finite di- 
mensional associative algebra. Set m = dim M(p). Let 
{pi, Pi, • • ■ j Pm} be a basis of the linear space M(p). 
Lemma 2. The metric tensor matrix fi, with entries 
given by tt(p)ij = Tr{p i p :j ), i,j = 1,2, ••• ,m, is non- 
singular. 

Proof. We prove by contradiction. If dct[f2(/?)] = 0, 
then the m row vectors of the matrix Q(p) are linear 
dependent. There exist ci, C2, • ■ • , c m G C which are not 
all zero, such that 

m 

^2 / c i Tr(p i p j )=Q, j = l,...,m. 
It follows that 

Tr((ci Pl +C2 P 2 H \-c m p m )Pj) = 0, j = l,...,m. (5) 

Since (A il A'j i ■ ■ ■ Ai r A'-^ G M(p) and pi are linear com- 
binations of the terms like A^A^ ■ ■ ■ Ai r Aj r , we have 
(cipi + c 2 p2 H h c m p m y G R{p) and 

(cipi + c 2 p2 H h c m p m ) t = hxpi + h 2 p 2 H h h mPm 

for some hi G C, i = 1, 2, ■ • • , m. It follows form (JS|) that 

Tr[(cipi+C 2 p 2 -\ \-C m Pm)(ciPl+C 2 p2^ l-CmPm^} = 0, 

which implies C\p\ + C2P2 + V c m Pm — 0. Hence one 

concludes that {pi, i = 1, 2, • ■ • , m} are linear dependent, 
which contradicts that {pi, p 2 , • • • , p m } is a basis of the 
linear space M(p). □ 

Actually, one can also prove that if {pi,i = 
1,2, ■•• , to} is a sequence of matrices such that 
det[fi(p)] ^ 0, then {pt,i = 1,2, ■•• , m} is linear inde- 
pendent. 

The invariants Tr[(AiAt)(A k Aj) ■ ■ ■ {A h Al)\ can be 

equivalently written as Tr[(A\Aj)(A\A{) ■ ■ ■ (A^Ap)]. 
Correspondingly one has the linear space defined by 

N(p) = spanHAlA^AlAt) ■ ■ ■ (A\A p )}, 

i,j,k,l,...,h,p — with finite dimension d = 

dimlN(p). If {0i, 2 ,-- - , d} is a basis of N(p) like 
R{p), we similarly have that the matrix with entries 
given by @(p)ij = Tr(0j0j), i,j — 1,2, ■•• ,d, satisfies 
det[6(»] ^ 0. 

Lemma 3. If two density matrices p = Y^ii=i ^i\ v i){ v i\ 
and p' = Y^i=i^i\ v 'i){ y 'i\ have the same values for the 
following invariants: 

Tr[{A l A]){A k A\)---{A h A% (6) 
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i,j,k,l,...,h,p = 1, ...,n, then the corresponding linear 
spaces M(p) and M(p') have the same dimension. More- 
over, the bases of M(p) and M(p') have a one to one 
correspondence: if (AiA^){A k A\) ■ ■ ■ (A^Aj,) is a basis el- 
ement of M(p), then (A' i A'J)(A' k A'^) ■ ■ ■ (A' h A'J) is a basis 
element of M(p'). Similar results hold between the linear 
spaces N(p) and M(p'). 

[Proof]. Assume that {p a , a — 1,2, ■ ■ ■ , m} is a basis of 
M{p). Each p a has a form of (A i A f j )(A k Al) ■ ■ ■ {A h Al). 
Denote p' a = {A'^A 1 ^ ) • ■ • (A' h A$), a = 
1,2, ■■■,m. Since det[Q(p)] = det[fi(p')] ^ from 
Lemma 1 and the condition ([6]), we have {p' a ,a = 
1,2, ■■ ■ , m} is linearly independent. Hence dim IR(p') > 
dim IR(p). In a similar way one can prove that 
dim iR(p) > dim M(p'). Therefore dim IR(p') = 
dimiR(p) and {p' a , a — 1,2, ■ ■ ■ , m} is a basis of M(p'). 

Similar results between the linear spaces IN(p) and 
IN(p') can be obtained from the expression of invariants 
TrUlA^AlAO ■ ■ ■ (A{A P )}. □ 

Local equivalence of bipartite states. — We now give 
the necessary and sufficient condition for the equivalence 
of bipartite states under local unitary transformations. 
Theorem Two arbitrary dimensional bipartite density 
matrices are equivalent under local unitary transforma- 
tions if and only if there exit eigenstate decompositions 
(H|) such that the following invariants have the same val- 
ues for both density matrices: 

.P{p)=Tr 2 (Tr lP s ), s = l,...,N 2 ; 
TrU^A^.-.iAnAl)], 

where i,j, k, I, h,p = 1, n. 

[Proof]. We have proved that the quantities in (JT)) are 
invariants under local unitary transformations. We now 
prove that these invariants are complete. Suppose that 
the states p = Y%=i and p' = Y%=i K\ u 'i)( v l\ 

have the same values to the invariants in ([7]). From 
J s (p) = J s (p'), s = 1,...,N 2 , we have that the two 
density matrices p and p' have the same set of eigenval- 
ues, i.e. \i = X'i for i = 1, 2, • • ■ ,n. 

Next we introduce the dual basis p* in M(p) such that 
Tr{ Pi p*) = Sij. In fact, let n(p)- 1 = [fi« (p)j, then 

m 

p* = J2^ ij iP)Pj- (8) 

3=1 

By Cramer's rule n(p)- 1 = \n(p)\- 1 Adj(ti)(p), so ft^(p) 
are given by polynomials in fiy (p) , consequently p* are 
expressed as polynomials in the basis element pi. 

Now we study the algebra structure of M(p). First we 
define the bilinear form ( , ) on M(p) by 

(<t,t) =Tr(ar+), a,r G M(p). (9) 

Then M(p) becomes a Hilbert space and each operator 
Pi is bounded under the norm ||p|| = (p, p} 1 / 2 . Next we 



have 

m 

PiPj = ^2 Tr (PiPjP*k)Pk- (10) 

k=l 

The invariance of the traces ([7]) implies that in the alge- 
bra M(p') the element p'* is given by the same formula 
as p*: 

m m 

P T=J2^(p')Pj=T, nlJ (P)Pr C 11 ) 

3=1 J = l 

Moreover p^p'^ are also given by the same structure con- 
stants: 

m m 

p'iP'i = £ Trip^p'M = £ Tr{p iPj pl)p' k . (12) 
fc=i fc=i 

Eqs (fTU)) and (fT2"|) imply that the map p^ — > p\ gives 
an isomorphism from the algebra -R(p) onto the alge- 
bra M(p'). Therefore M(p) and M(p') are two equivalent 
representations of the same underlying associative alge- 
bra. Thus there exists a non-singular matrix T such that 
pi = Tp^T^ 1 for all i = 1, . . . , m. In particular, we have 
A. L A\ = TA'iA'^T- 1 , i = 1, . . . ,m. As AiA\ are hermi- 
tian, due to the algebraic property ((§]) and using Theorem 
12.36 in [H we have 

AiA^uA'^u^ (13) 

where u is the unitary part of T in the polar decomposi- 
tion. 

Similarly from isomorphism of the algebras iV(p) ~ 
IN(p'), one has 

A\A i =wA!}A! i w\ (14) 

for some unitary matrix w and all i. 

By using relations ([IB")) and (TT4"|) we can show 

A, =it^O*)*, i= l,...,ri. 

In fact, let Ui and be the unitary matrices that 
diagonize the hermitian matrices AiA\ and A\A.- L re- 
spectively, UiA t A\u\ = diagir] 2 ^...,^}, u^A^ = 
diaglrj^, ...,rjf }. From the procedure of the singular- 
value decomposition of matrices j2lj |. we have UiAiU 1 ^ — 
diagirj^, ...,i]i N } with r\ k > 0. From $13§ and (fl4| . 
we have UjAjAjuj = u^A^A'^u^u] = diagln^ , rft N } 
and u^AjA^ 1 = tt>A< + A>V.t = ^{77? , ...,»&,}. 
Therefore from singular value decomposition we have 
UiuA^w^u'^ = diag{j]i x , ...,rji N }. Hence we obtain 
^Aiu'^ = UiuA'^u'^, i.e. A l = uA'^ = uA'^w*) 1 . 

Since A l = uA'^w*) 1 , we have = ® (w*)^^), 
i = 1, n, and p' = ®(w*)^ p u<E>w* . Therefore p' and 
p are equivalent under local unitary transformations. □ 
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From the proof one sees that the invariants (J7J are 
complete, finite and can be easily calculated. For non- 
degenerate states p and p', the eigenvector of p cor- 
responds uniquely to the eigenvector |z/|) of p' . Hence 
the correspondence between A4 and A[ is also unique. It 
is straightforward to judge the local unitary equivalence 
of p and p' by simply comparing the values of invariants 
from Ai and A\ one by one. 

For degenerate states, it becomes less operational since 
the eigenvector decompositions ((T|) is no longer unique. If 
\vi), i = 1, r, are the orthogonal eigenvectors with re- 
spect to a same eigenvalue, then \vA = J2j ^ij\ v j) 1S a ^ so 
a set of linearly independent orthogonal eigenvectors to 
the same eigenvalue for any unitary matrix (U)ij = Uij. 
Therefore the expressions of the invariants are not 
unique in general. However, by linear combinations of 
the invariants in ([7]) we can get a set of invariants which 
does not depend on the detailed eigenvector decomposi- 
tions, 



^2i,i',j, 3 ',...,k,k'=i Tr{AiA\ t AjA^, ■ 
Yji,i<,j,j<,...,k,k'=i Tr(A i A i ,A ;j Aji ■ 



■A k Al), 



A\Ak. 



(15) 



where the indices k'} is any given permutation 

of the indices k}. Different given permutation 

gives different invariants. Since the invariants (|15j) are 
independent on the detailed eigenvector decompositions, 
they gives an operational necessary condition for the lo- 
cal equivalence, with respect to the subspaces spanned by 
the degenerate eigenvectors. If two density matrices with 
the same degenerate eigenvalue are equivalent under local 
unitary transformations, they must have the same values 
of all the invariants in (|15|) . As an simple example, con- 
sider p = diag{l/2, 1/2,0,0}, p' = diag{l/2, 0, 1/2, 0}. 



We have A x = A[ = 



1 




, a 2 = (A> 2 y = 



1 





From (|15j) we have that p and p' are not local uni- 
tary equivalent, since Tr(AiA\A x A\ + A x A\A 2 A\ + 
A 2 A\A 1 A\ + A 2 A\A 2 A\) = 2, while Tt{A\A^_A\A^ + 
A 1 A 1 4't j_ A' 4't A 1 -l- A' A' A'^\ — A 

A \ A 2 A 2 A \ ' /i 2 /1 l A \ A 2 + A 2 A 2 A 2 A 2 ) ~ 

Conclusion and remarks. - We have investigated 
the nonlocal properties of arbitrary dimensional bipartite 
quantum systems and solved the local equivalence prob- 
lem by presenting a complete set of invariants such that 
two density matrices are locally equivalent if and only 
if all these invariants have the equal values. Although 
the independent invariants may vary with the detailed 
bipartite states, the number of the invariants one needs 
to check is finite. Here we have dealt with the case that 
the dimensions of both Hilbert spaces are the same for 
simplicity. Nevertheless the case that the Hilbert spaces 
have different dimensions can be similarly discussed. Our 
approach of constructing local invariants and algebraic 
proof of the sufficiency may also shed light on multipar- 
tite case, for which only the multipartite pure state case 
has been extensively studied 3, ljjj . 
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